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O ' Abstract 

(n: 

Transport of a particle in a spatially periodic harmonic potential under the influence of a slowly 

W ', time-dependent unbiased periodic external force is studied. The equations of motion are the same as 

^ ' in the problem of a slowly forced nonlinear pendulum. Using methods of the adiabatic perturbation 

•i-H ' 

theory we show that for a periodic external force of general kind the system demonstrates directed 
(ratchet) transport in the chaotic domain on very long time intervals and obtain a formula for the 
average velocity of this transport. Two cases are studied: the case of the external force of small 

OS : 

amplitude, and the case of the external force with amplitude of order one. The obtained formulas 
^. . can also be used in case of a non-harmonic periodic potential. 
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I. INTRODUCTION 

In recent years, studies of transport phenomena in nonlinear systems have been attracting 

a growing interest. In particular, a large and constantly growing number of papers are 

devoted to dynamics in systems which allow for directed [on average] motion under unbiased 

external forces and are referred to as ratchet systems. Their intensive study was motivated 

by problems of motion of Brownian particles in spatially periodic potentials, unidirectional 

transport of molecular motors in biological systems, and recognition of "ratchet effects" 

in quantum physics (see review |l| and references therein). Generally speaking, ratchet 

phenomena occur due to lack of symmetry in the spatially periodic potential and/or the 

external forcing. It is interesting, however, to understand microscopic mechanisms leading 

to these phenomena. A possible approach is to neglect dissipation and noise terms and 

arrive at a Hamiltonian system with deterministic forcing. Thus, one can make use of 

results obtained and methods developed in the theory of Hamiltonian chaos. Many papers 

studying chaotic transport in such Hamiltonian ratchets appeared in the last years (see, 

e.g. [2|, y, U, la, la, lit Iq])- In particular, in [6| the ratchet current is estimated in the case 

when there are stability islands in the chaotic domain in the phase space of the system. The 

(1 
borders of such islands are "sticky" [9( and this stickiness together with desymmetrization 

of the islands are responsible for occurrence of the ratchet transport. 

We consider the problem of motion of a particle in a periodic harmonic potential U(q) = 
UQCOsq, where q is the coordinate and u = const, under the influence of unbiased time 
periodic external forcing. The equations are the same as in the paradigmatic model of a 
nonlinear pendulum under the action of external torque with zero time average. We study 
the case when the external forcing is a time-periodic function of large period of order e -1 , 
< £ -C 1, and use results and methods of the adiabatic perturbation theory. If e is small 
enough, there are no stability islands in the domain of chaotic dynamics (see [10 ). Thus, 
the mechanism of ratchet transport in this system differs from one suggested in |6[. 

The main objective of this paper is to find a formula for the average velocity V q = (q) 
of a particle in the chaotic domain on very large time intervals. Chaos in the system is a 
result of multiple separatrix crossings due to slow variation of the external forcing. At each 
crossing the adiabatic invariant ("action" of the system) undergoes a quasi-random jump 
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121] ). We show that these jumps result in effective mixing and uniform distribution 



of the action along a trajectory in the chaotic domain. On the other hand, direction and 
value of velocity depends on the immediate value of the action. Thus, to find the average 
velocity of transport on time intervals of order or larger than the mixing time, we find 
formulas for displacement in q at a given value of the action and then integrate them over 
the interval of values of the action corresponding to the chaotic domain. We demonstrate 
that for an external force of general kind (i.e. with zero time average but lowered time 
symmetry, cf. [5|), there is directed transport in the system and obtain an analytic formula 
for the average velocity V q of this transport. 

In Section 2, we obtain the main equations in the case of external force of small amplitude 
and describe the diffusion of adiabatic invariant due to multiple separatrix crossings. This 
diffusion makes dynamics chaotic in a large domain of width of order 1. In Section 3, we 
derive the formula for the average velocity of transport on very long time intervals (of order 
of typical diffusion time and larger) and check it numerically. In Section 4, we consider the 
case when the external forcing is not small, of order 1. Width of the chaotic domain in 
this situation is large, of order e~ 1 . In this case chaos arises as a result of scattering of the 
adiabatic invariant on the resonance. Formula for the average velocity of transport in this 
case turns out to be much simpler than in the case of external forcing of small amplitude. 

II. MAIN EQUATIONS. DIFFUSION OF THE ADIABATIC INVARIANT 

We start with a basic equation describing a nonlinear pendulum under the influence of 
an external force f(t): 

q + u%sinq = f(t). (1) 

We assume that /(£) is small in amplitude periodic function of time of large period: f(t) = 
ef(r) = ef(r + 2n), where r = et. Moreover, we assume that f(t) is a function with zero 
time average. This system is Hamiltonian, with time- dependent Hamiltonian function 

2 

H = ^--u 2 cosq-ef(et)q. (2) 

We make a canonical transformation of variables (p, q) h-» (p, q) using generating function 
W = (p—F(et))q, where F(r) = — f /(r)dr. Thus, F(r) is a periodic function defined up to 
an additive constant, which we are free to choose. To make the following presentation more 
clear, we choose this constant in such a way that the minimal value of F is F min > Aloq/iz. 



Note that q = q. After this transformation, Hamiltonian of the system acquires the form 
(bars over q are omitted): 

in n-rW-i 

(3) 



(p-F(r)f 2 
ii = Uq cos g. 



Phase portrait of the system at a frozen value of r (we call it the unperturbed system) is 
shown in Fig. [TJ The separatrix divides the phase space into the domains of direct rotations 
(above the upper branch of the separatrix), oscillations (between the separatrix branches), 
and reverse rotations (below the lower branch of the separatrix). Introduce the "action" 
I associated with a phase trajectory of the unperturbed system on this portrait. In the 
domains of rotation, I equals an area between the trajectory, the lines q = —n, q = tt, and 
the axis p = 0, divided by 2tt; in the domain of oscillations, this is an area surrounded by 
the trajectory divided by 27r. It is known that I is an adiabatic invariant of Q: far from 
the separatrix its value is preserved along a phase trajectory with the accuracy of order e 
on long time intervals (see, e.g., }13|). 




FIG. 1: Phase trajectories of the system ([3]) with sample function F(r) = 
A (l + 2exp [— a(sinr) 2 ]), A = 10, a = 8, r = ir/2, loo = 1. The bold line is the separa- 
trix. 

Location of the separatrix on the (q, p)-plane depends on the value of F(t). As t slowly 
varies, the separatrix slowly moves up and down, and phase points cross the separatrix 
and switch its regime of motion from direct rotations to reverse rotations and vice versa. 
Recall known results on variation of the adiabatic invariant when a phase point crosses the 
separatrix. The area surrounded by the separatrix is constant, and hence, capture into the 
domain of oscillations is impossible in the first approximation (in the exact system, only a 
small measure of initial conditions correspond to phase trajectories that spend significant 



time in this domain; thus their influence on the transport is small). To be definite, consider 
the situation when the separatrix on the phase portrait slowly moves down. Thus, phase 
points cross the separatrix and change the regime of motion from reverse rotation to direct 
rotation. Let the action before the separatrix crossing at a distance of order 1 from the 
separatrix be I = J_ and let the action after the crossing (also at a distance of order 1 
from the separatrix) be I = I + . In the first approximation, we have J + = /_ + 8uJo/n, i.e. 
the action increases by the value of the area enclosed by the separatrix divided by 2tt (see, 
e.g., [liJ, ll5|). We shall call this change in the action a "geometric jump". If the separatrix 
contour slowly moves up, and a phase point goes from the regime of direct rotation to the 
regime of reverse rotation, the corresponding value of the action decreases by the same value 
8u)o/tt. Thus, in this approximation, the picture of motion looks as follows. While a phase 
point is in the domain of reverse rotation, the value of I along its trajectory stays constant: 
I = I_. After transition to the domain of direct rotation, this value changes by the value 
of the geometric jump. The transition itself in this approximation occurs instantaneously. 
After the next separatrix crossing, the adiabatic invariant changes again by the value of the 
geometric jump, with the opposite sign, and returns to its initial value /_. We call this 
approximation adiabatic. 

In the next approximation, the value of action at the separatrix crossing undergoes a 
small additional jump. Consider for definiteness the case when the separatrix contour on 
the phase portrait moves down, and J_ and I + are measured when it is in its uppermost 



and lowermost positions, accordingly. Results of III, [12] imply the following formula for the 
jump in the adiabatic invariant: 

27r(J+-J_) = 16w + 2o(l - £)e01n(ee) 

+ agein 27r ^"^ -26ee(l-0, 

r (0 

where a = Uq 1 , b = UqHii^uJq) , 6 = 2tcF'(t*). Here F' is the r-derivative of F, T* is 
the value of r at the separatrix crossing found in the adiabatic approximation, T(-) is the 
gamma-function. Value £ is a so-called pseudo-phase of the separatrix crossing; it strongly 
depends on the initial conditions and can be considered as a random variable uniformly 
distributed on interval (0, 1) (see, e.g., [12j). Thus, value of the jump in the adiabatic 
invariant at the separatrix crossings has a quasi-random component of order sine. 

Accumulation of small quasi-random jumps due to multiple separatrix crossings produces 



diffusion of adiabatic invariant (see, e.g., 13]). On a period of F(r) (after two separatrix 
crossings) the action changes by a value of order elne. Hence, after iV ~ e~ 2 (lne)~ 2 separatrix 
crossings the adiabatic invariant varies by a value of order one. As a result, in time of order 
e~ 3 (lne)~ 2 the value of adiabatic invariant is distributed in all the range of values 



tdiff 

corresponding to the domain where phase points cross the separatrix on the phase plane; its 
distribution is close to the uniform one. We have checked this fact numerically for a sample 
function F(t) at various parameter values. Poincare sections and distribution histograms of 
I in all the cases look similar; see an example in Fig. [2j 
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FIG. 2: Left panel: Poincare section at r = 0mod27r of a long phase trajectory (5 • 10 4 dots) . 
All the points are mapped onto the interval q £ (— ir, n). F(t) = A (l + 2exp [— a(sinT) 2 ]) with 
A = 10, a = 16, £ = 0.005, uj$ = 1. The empty region in the chaotic sea corresponds to phase 
points eternally locked in the domain of oscillations; they never enter the chaotic domain and do 
not participate in the transport. Right panel: Histogram of / on the segment (I m in, Imax —8u>o/ir)) 
along the same phase trajectory. 



III. AVERAGE VELOCITY OF THE TRANSPORT 



Our aim is to find a formula for average velocity V q along a phase trajectory on time 
intervals of order tdiff or larger. We first only take into consideration the geometric jumps, 
and afterwards, to obtain the final result, we take into account the mixing due to small 
quasi-random jumps. To simplify the consideration, assume that function F(t) has one 
local minimum F min and one local maximum F max on the interval (0, 2ir). The main results 
are valid without this assumption. 



Introduce /, denned in the domains of rotation as follows: it equals the area bordered 
by the trajectory, the line p = F(t), and the lines q = —n, q = tt, divided by 2n. Thus, 
I = \F(t) — I\. Frequency of motion in the domains of rotation is u(I), where u>(I) at 
I > Auq/tt is the frequency of rotation of a standard nonlinear pendulum with Hamiltonian 
Hq = p 2 /2 — oJq cosq, expressed in terms of its action variable /. We do not need an explicit 
expression for function u(I). From Hamiltonian ([3]) we find q = p — F(t). Consider a phase 
trajectory of the system frozen at r = f in a domain of rotation. Let the value of action on 
this trajectory be I = I Q . Then the value of q averaged over a period T of rotation equals 

fi\q\dt/T = 2n/T = uj(\F(T)-Io\). 

Now consider a long phase trajectory in the case of slowly varying r. Let on the interval 
(ti, t 2 ) a phase point of ([3]) be below the separatrix contour. In the adiabatic approximation, 
the value Iq of the adiabatic invariant along its trajectory is preserved on this interval. Hence, 
at r G (7*1, t 2 ) we have 

2vrF(r) - 2ttI > 8u , (4) 

and the equality here takes place at r = T\ and t = T2. In the process of motion on this 
time interval, q changes (in the main approximation) by a value 

Ag_(/ ) = -- r l u(F(r) - J )dr. (5) 

On the interval (r 2 , T\ + 2ir) the phase trajectory is above the separatrix contour, and the 
value of the adiabatic invariant equals 1$ = Iq + 8cjo/vt due to the geometric jump. On this 
interval we have 

2vrF(r) - 2ttI < 8u . (6) 



In the process of motion on this time interval, q changes by a value 

Ag + (/ ) = - [" "^ u(\F(r) - I \)dr. (7) 



^ * n 



Total displacement in q on the interval (r l7 Ti +2n) equals Ag(J ) = Ag_(J ) + Aq + (I ), and 
the average velocity on this interval is eAq(Io)/(2Tr). 

Consider now the motion on a long enough time period At ~ tdiff- Due to the diffusion in 
the adiabatic invariant described above, in this time period values of Iq, defined as a value of 
/ when the phase point is below the separatrix contour, cover the interval (/ m » n , I max —8uo/7i). 
Here I min = F min — 4u /ti and I max = F max + 4u /n. Assuming that the distribution of / 



on this interval is uniform, to find the average velocity, we integrate eAq(Io)/(27r) over this 
interval. Integrating (jSJ) over Iq and changing the order of integration we find 



flmax c 


ujq/tt 


-i p27T /•F(t)—4luo/it 


/ 


Ag_d/ =- 


— / dr / tu(F(r) - I )dl 


J-min 




£ Jo Ji min 

i p2tt rF{T)—I min 




= 


— / dr / Lu(rj)drj. 

£ JO J^o/lT 



p2tt rF(r)-I mirl r2ir 

- dr u( V )dr, = - (H (F(r)-I mi/n )-H s )dr, (8) 

Jo JaujcI-k Jo 



Now we take into account the equality to(I) = 8Hq(I) / dl (recall that Hq(I) is the Hamil- 
tonian of a nonlinear pendulum as a function of its action variable) and obtain 

r°2n rF(r)-I mirl r2ir 

dr / Lu(r])drj = — 

Jaluo/tt Jo 

where Hq is the value of H on the separatrix. Similarly, integrating (|7j) we obtain 

Adding (jHJ) to (j2j) and dividing by 27i(F max — F min )/e we find the expression for the average 
velocity V q of transport on long time intervals: 

1 
q ~ 2it(F - F ■ N 

^ n y- 1 - max ± mi 



f 



In i p2ir 

Aq + dl = - (H (I max -F(T))-H°)dT. (9) 

£ Jo 



mm 

2tt 



i>Ztt 

x / (H (I max -F(r))-H (F(r)-I mm ))dT. (10) 

In ( TTU1) . Hq(I) can be found as the inverse function to 1(h), which defines action as a 
function of energy in domains of rotation of a nonlinear pendulum. For the latter function, 
the following formula holds (see, e.g., 161]): 

4 
1(h) = -u q kE (l/«), k> 1, (11) 

7T 

where k 2 = (l + h/u> i )/2, £(■) is the complete elliptic integral of the second kind. If function 
F(r) has several local extremes on the interval (0, 2ir), F min and F max in ffTUj) are the smallest 
and largest values of F respectively. 

It can be seen from ffTUl) that for function F(t) of general type V q is not zero, and hence 
there is the directed transport in the system. We checked this formula numerically for a 
sample function F(r) = A (1 + 2exp [— a(sinr) 2 ]) , a > at various values of parameters e 
and a. Typical plots of q against time t are shown in Fig. [31 



8 





a = 4 


a = 8 


a = 16 


£ = 0.1 


4.721 


6.756 


8.363 


£ = 0.05 


4.446 


6.681 


8.076 


£ = 0.01 


4.298 


6.211 


7.442 


£ = 0.005 


4.598 


6.702 


8.202 


jrtheor 


4.393 


6.679 


8.110 



TABLE I: Numerically found values of V q corresponding to various values of parameters e, a (four 
upper rows, A = 10, loq = 1) and theoretical values Vg heor obtained according to (|10p (the bottom 
row). 
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FIG. 3: Left panel: q against t for ten different initial conditions (comparatively short time interval), 
a = 4, £ = 0.05. Central panel: q against the number of periods of the external force for a sample 
trajectory (10 4 periods), a = 16, £ = 0.05. Right panel: q against the number of periods of the 
external force for the same trajectory (10 5 periods), a = 16, e = 0.05. Parameter A = 10 in all 



cases. 



The results of numerical checks of formula (fTUj) are represented in Table [H To find 
numerical values of V q presented in the table, we integrated the system with Hamiltonian 
([3]) on a long time interval At = 2rr ■ 10 6 /£: with a constant time step of 7r/100 (5th order 
symplectic scheme [17]). Use of a symplectic scheme for long time simulations of Hamiltonian 
systems is necessary in order to ensure that creeping numerical error do not end up washing 
off the invariant tori bounding the chaotic domain. The table demonstrates satisfactory 
agreement between the formula and the numerics. 

Finally, we note that formula ( fTUl) can be used also in the case of arbitrary (non-harmonic) 
spatially-periodic time independent potential in place of the term — cugcosg in (J2J), (J3J). Of 



course, in this case function Hq is different from the Hamiltonian of the nonlinear pendulum, 
but it always can be found, at least numerically. 

IV. THE CASE OF NOT SMALL EXTERNAL FORCING 

In this section we study the case when the external forcing is not small. In this case, 
amplitude of function f(t) in ([T]) is a value of order one: /(£) = /(r) = f(r + 2n). The 
equations of motion are: 

Q = P, P = -^o sin ^ + /( r )' f = e - ( 12 ) 

One can see from the second equation, that magnitude of momentum p can reach values 
of order e -1 . Making the canonical transformation with generating function W = (p — 
e~ 1 F(et))q, where F(t) is defined in Section II, we obtain the Hamiltonian: 

H= {p - £ ~ lF{T))2 -^cos q . (13) 

Introduce p = ep and rescaled time t = e~ x t. We denote the derivative w.r.t. t with prime 
and thus obtain: 

q' —p — F(t), p = -e 2 u 2 sin q, t' = e 2 . (14) 

This is a system in a typical form for application of the averaging method. We average over 
fast variable q and obtain the averaged system 

p = i T ! = e 2 . (15) 

The averaged system describes the dynamics adequately everywhere in the phase space 
except for a small neighborhood of the resonance at p — F(t) = 0, where the "fast" variable 
q is not fast. When a phase trajectory of the averaged system crosses the resonance, value 
of the adiabatic invariant (in this case it coincides with p) undergoes a quasi-random jump 



of typical order ye^ = e (see, e.g., 13|, |l8j). Thus, the situation in this case is similar to 
one studied in Section 2, but the typical value of a jump is of order e, and, accordingly, the 
typical diffusion time is tdiff ~ £~ 3 - Phase trajectories of the averaged system that cross the 
resonance correspond to values of p belonging to the interval (F m i n , F max ). Therefore the 
chaotic domain of the exact system is, in the main approximation, a strip F m j„ < p < F max . 
Captures into the resonance followed by escapes from the resonance (see [13j, ll8() are also 
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possible in this system. However, probability of capture is small, of order e, and hence 
impact of these phenomena on the transport is small. 

In Fig. H]we represent a sample of Poincare section of a long phase trajectory of (1131) and 
the corresponding histogram of p for this trajectory. The plots show that the distribution 
of p is close to the uniform one. 
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FIG. 4: Left panel: Poincare section at r = 0mod27r of a long phase trajectory (5 • 10 4 
dots) of system (|13p . All the points are mapped onto the interval q S (— tv, 7r). F(t) = 
A (l + 2exp [— o(sinr) 2 ]) with A = 0.5, a = 4, e = 0.01, ujq = 1. Right panel: Histogram of p 
along the same phase trajectory. 



Like in Section 3, one can find value Aq on one period of perturbation, then average over 
the range of adiabatic invariant corresponding to the chaotic domain, and find the average 
velocity of transport in this case. Thus we find: 



Aq 



2tt/£ 



pdt 



" v^u ! 



t ./0 



2~ 



(p-F(r))dr. 



(16) 



To find V q , we have to integrate this expression over p from F min to F max (i.e., over the 
chaotic domain) and divide the result by (F max — F min ) and by the length of the period of 
the external forcing 2^/e. Thus we obtain 

* £* max 



Vn 



e 



Dm i -L'min) J F m i n 



Aqdp. 



q 2tt(F — F 
Substituting Aq from f[T6"j) and integrating, one straightforwardly obtains 



(17) 



V q = ^- I (F max + F rnm -2F(r))dr. 



(18) 
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a = 1 


a = 2 


a = 4 


£ = 0.1 


0.046 


0.128 


0.253 


e = 0.05 


0.046 


0.112 


0.225 


£ = 0.01 


0.0353 


0.1050 


0.2044 


e = 0.005 


0.0369 


0.1081 


0.1916 


-irtheor 

t V q 


0.0389 


0.1018 


0.2006 



TABLE II: Numerically found values of eV q corresponding to various values of parameters e, a 
in system (fl~3j) for F(r) = A (1 + 2exp [— a(sinr) 2 ]) (four upper rows, A = 0.5, ujq = 1). In the 
bottom row theoretical values eV^ heor obtained according to (|19p are shown. 



Another possibility to find V q is to use already obtained formula ffTOj) . This way leads to the 
same result. When using (fTUj) . one should keep in mind that in the considered case H in 
this formula is the Hamiltonian of the averaged system, i.e. H Q = I 2 /2, I max = F max , I min = 
F m i n , and that according to (Tl~3l) we should put e~ l F instead of F everywhere in the formula. 
Thus, ( TTUl) is much simplified, and we again arrive at formula ( TT8l) . Factor e" x is due to the 
fact that a phase point spends significant time moving at large velocities corresponding to 
p ~ e^ 1 . 

Note that formula ffl~8]) can be rewritten in a more elegant form as: 

J- / ± max T J- mi; 



V n 



(F(r)) 



(19) 



s v 2 

where the angle brackets denote time average. The results of numerical checks of the formula 
are represented in Table [Til Like in the previous Section, we integrated the system with 
Hamiltonian (Tl~3l) on a long time interval At = 2tt ■ 10 6 /£. 

Remarkably, formula f|T9|) is valid for arbitrary smooth periodic potential (not necessar- 
ily harmonic) in place of the term — Uq cos q in ([2]) , ([3]) . The potential may also depend 
periodically on time with the same period as that of the external force. 



V. SUMMARY 

To summarize, we have considered the phenomenon of the directed transport in a spatially 
periodic harmonic potential adiabatically influenced by a periodic in time unbiased external 
force. We have shown that for the external force of a general kind the system exhibits 
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directed transport on long time intervals. Direction and average velocity of the transport 
in the chaotic domain are independent of initial conditions and determined by properties 
of the external force. We studied two different cases: the case of small amplitude of the 
external force and the case, when this amplitude is a value of order one. We have obtained 
an approximate formula for average velocity of the transport and checked it numerically. 
The final formulas (ITU]) and (TT9I) are valid for any smooth periodic potential (not necessarily 
harmonic one). 
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